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We analyze the dilepton and photon emission rates from a 
hadronic gas using chiral reduction formulas and a virial ex- 
pansion. The emission rates are reduced to pertinent vacuum 
correlation functions, most of which can be assessed from ex- 
periment. Our results indicate that in the low mass region, 
the dilepton and photon rates are enhanced compared to most 
of the calculations using chiral Lagrangians. The enhance- 
ment is further increased through a finite pion chemical po- 
tential. An estimate of the emission rates is also made using 
Haag's expansion for the electromagnetic current. The rele- 
vance of these results to dilepton and photon emission rates 
in heavy-ion collisions is discussed. 



1. Recent fixed target experiments at the CERN 
SPS using ultra-relativistic heavy-ion collisions have re- 
ported an excess of dielectrons (CERES) Q and dimuons 
(HELIOS3) j| over a broad range of dilepton invariant 
mass starting from the two-pion threshold. Some en- 
hancement in the photon emission rate has been reported 
by WA80 H although the direct photon measurements 
by CERES and HELIOS3 seem to suggest that the ex- 
cess is within statistical errors when the Dalitz decays are 
subtracted out. The latter constitute 90% of the photon 
spectrum below the two-pion threshold. 

Enhancements in either the dilepton or photon spec- 
tra do not seem to follow from calculations using chiral 
Lagrangians in hadronic phase [ffl|5) or high temperature 
QCD for the quark-gluon phase ||, although neither of 
them is assumption free for temperatures around the pion 
mass. A rate departure from p-A collisions may be in- 
dicative of some medium modifications in the hadronic 
phase, although the rate may be enhanced or depleted ||] 
depending on model assumptions. 

In view of this, it is important that the rate calcu- 
lations are reassessed in as model independent way as 
possible. Also, since the dilepton and photon emission 
rates originate from the same current-current correlation 
albeit with different kinematics, it is important that both 
rates are assessed simultaneously in a consistent man- 
ner. Better photon and dilepton measurements are yet 
to come, putting to test our understanding of the strong 
interactions in the hot hadronic phase. 

In this letter we will show that in a thermal hadronic 
environment with a temperature near or below the pion 
mass the dilepton and photon emission rates are directly 
related to the absorptive parts of the scattering ampli- 
tudes of virtual and real photons off an arbitrary number 



of on-shell pions. These scattering amplitudes are re- 
ducible to pertinent vacuum correlation functions using 
chiral reduction formulas || . To leading order in the pion 
density, the emission rates are constrained by data from 
electro-production and -annihilation, tau-decays, pion ra- 
diative decay and two-photon fusion reactions. Our re- 
sults indicate that in the low mass region the dilepton 
and photon rates are enhanced in comparison with cal- 
culations using explicit reaction processes from chiral La- 
grangians. There is also an enhancement for finite pion 
chemical potential in both cases. For comparison, we 
also give a direct estimate for the emission rates using 
Haag's expansion for the electromagnetic current. The 
relevance of our results to lepton and photon emission 
from heavy-ion collisions is discussed. 

2. In a hadronic gas in thermal equilibrium, the rate 
R of dileptons produced in a unit four volume follows 
from the thermal expectation value of the electromag- 
netic current-current correlation function JTcj ] . For mass- 
less leptons with momenta p\,pi, the rate per unit dilep- 
ton momentum q = pi + P2 is given by 
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where a = e 2 /4ir is the fine structure constant, 
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M is the hadronic part of the electromagnetic current, 
H is the hadronic Hamiltonian, F the free energy, T 
the temperature, and the trace is over a complete set of 
hadron states. We have set the pion chemical potential 
to zero. The results for a non-zero pion chemical poten- 
tial are discussed in section 5. For leptons with mass mi, 
the right-hand side of (Q) is multiplied by 
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From the spectral representation and symmetry, the 
rate may also be expressed in terms of the absorptive 
part of the time-ordered function 



W(q) = 
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From here on we take T < m, and consider only pion 
states. Expanding the trace in terms of incoming pions, 
and summing over disconnected pieces leads to the virial 
expansion 

W F (q)=i J d A xe iq - x <0|T*J' i (a;)J Al (0)|0> 

x < 7T a (fc) in| T*J^(x)J Al (0)| 7 r a (fc) in > colm . 

,J_V^ f ^ 3fc i d 3 k 2 n(uj kl ) n(uk 2 ) 
91 2^ 
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with the pion energy uj k = y/ m 2 + k 2 , and the ther- 
mal pion distribution n(u>) = l/(e"/ T — 1). The ma- 
trix elements in (|^) correspond to the forward scatter- 
ing amplitudes of a virtual photon off on-shell pions, 
7* wk —> 7* wr. 

To leading order in the pion density, 

d 3 k . . 

Wf n{uJk) > 

only the first two terms in (||) contribute. A slight re- 
finement is to note that the reduction of one pion is as- 
sociated with a factor of 1//^, where f v — 93 MeV is the 
pion decay constant, so that the dimensionless expansion 
parameter should be k = n/2m 7r f 2 . According to Ref. 
|Ll| , a pion gas with a temperature T < is character- 
ized by a pion mean free path of the order of 2-3 fm over 
a wide range of momenta, with a pion mean distance of 
the order of 2 fm. This gives k = 0.3, so the truncation 
should be reasonable. 

The absorptive part of the vacuum contribution is di- 
rectly related to data through e + e _ annihilation into 
hadrons The forward scattering amplitude 7*7r — > 
7*7r is unfortunately not measurable. However, it can be 
constrained by data as we now discuss. 

3. Decomposing the electromagnetic current into the 
isovector part V 3 and the isoscalar part B, the forward 
7*7r — > 7*7r matrix element involves three types of time- 
ordered correlators BB, BV, and VV. In pionic states, 
the BB and BV correlators are expected to be small. 
Indeed, in the soft pion limit, conventional PCAC gives 
zero for both matrix elements. The former because B 
commutes with the axial charge, and the latter because 
the vacuum is isospin invariant. The VV correlator in 
the one-pion state W F can be reduced to pertinent vac- 
uum correlators using chiral reduction formulas || . The 
result is (u = k ■ q) 
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+ {q^ -q) + {k -> —k) + (q, k -> -q, -k) (7) 
where A R (k) is the retarded pion propagator 
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j A is the axial vector current with the one-pion contri- 
bution subtracted, V is the isovector current discussed 
above, and a is a scalar density related to qq 13] . 11^ and 
Hv are the transverse part of the axial < 0|T*j^j^i|0 > 
and vector < 0|T*VV|0 > correlators respectively. The 
absorptive part of Ily follows from electroproduction 
data, while the absorptive part of 11^ follows from tau 
decay data. 

Note that the longitudinal part of the axial correlator 
does not contribute here. Taking the chiral limit (m, = 
0) and for a zero-momentum pion, 

Wf (g, 0) = 12g 2 Im (U v (q 2 ) ~ U A (q 2 )) 



which reproduces the results obtained in |lj,|13|]. As al- 
ready indicated, eqs. only takes care of the isovec- 
tor part of the electromagnetic current, both at zero and 
first order in the virial expansion. The inclusion of the 
isoscalar contribution to the zeroth order part leaves our 
arguments unchanged (at the 1% level) except around 
the phi region. 

The third term of (Q) contains the principle value in- 
tegral from (||). However, there is no pole within the 
range of integration unless q 2 = 2m 7r qo. This term is 
proportional to q 2 lmllv (q 2 ) vanishing at q 2 = 0. 

The last three terms of (|7|) are vacuum correlators 
which appear in other low-energy processes and are also 
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constrained by data. The j^j^V correlator appears in 
the radiative decay of the pion |[J . Empirically, the con- 
tribution of this term is small at threshold. There are 
two ways in which this can be checked analytically: a 
resonance saturation JL4| and a one- loop expansion f|. 
Both of these methods give zero for this correlator. 

The j/UviVV and VV<r correlators appear in the two- 
photon fusion process || . Resonance saturation for both 
of these gives on the order of a 1% correction and the 
one-loop result contributes about 1% in the low mass 
region and even less at higher energies. Both descriptions 
are in agreement with data from 77 — » 7r + 7r~ Jlq] and 
77 — > 7T 7r° Jig ]. Hence, the last three correlators in (Q) 
may be ignored and the rest can be extracted from data. 
This information will be used in section 5 to assess the 
dilepton and photon rate. 

The present analysis can be readily carried out to the 
photon emission rate through 



q d 3 q 4tt 2 W 
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with q 2 = 0. The constraints of current conservation 
and chiral symmetry apply equally well to the photon 
emission rate. As indicated above, it is important to 
assess both emission rates simultaneously for consistency. 
In our case, the only non-zero contribution comes from 
the II a term in eq. ([?]). 

4. For completeness, we mention another expansion 
scheme called the Haag expansion. From the pion elec- 
tromagnetic form factor 

< 7r*(pV(p) in|SJ A1 (x)|0> = 

= ie a3b (p' -p) p F v ((p + p') 2 ) 

and the crossed versions, we obtain the Haag expansion 
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and similarly for J M (y)S^, where S is the hadronic 
S-matrix. Inserting this into (Q), and using energy- 
momentum conservation along with Wick's theorem, we 
obtain 
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The form factor Fy(g 2 ) for g 2 > 4m 2 is directly mea- 
sured from electro-production. For q 2 < it follows from 
pion scattering on hydrogen targets. In both cases, rho 
dominance holds to a good accuracy. This is similar to 
the first term in the virial expansion (|), but the Haag 
expansion is presumably less controlled. In particular, 
the photon rate is zero in this approximation. We note 
that ([h]) in the q = limit has already been used in [H . 
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FIG. 1. The total integrated rate for dielectrons from a 
hadronic gas of T = 50, 100 and 150 MeV. The solid line is 
from the virial expansion and the dotted line from the Haag 
expansion. 



5. In Fig. 1, we show the numerical results follow- 
ing from the above analysis for the dielectron rate as a 
function of the invariant dielectron mass M — ^/q 2 , all 
the way up to the al mass for the temperatures T = 50, 
100, 150 MeV. The change of variables to the rapidity y 
and magnitude of the transverse momentum q± jl(| 
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and an integration over y and q ± is performed. The solid 
curve is the result following from (Q) using (0) and <^). 
The dotted curve is the result following from ([l]) and the 
Haag expansion (|l(]). For comparison, we also show the 
integrated dimuon rate in Fig. 2. As the temperature 
increases the rate for dilepton momentum below the rho 
pole increases with respect to the Haag expansion. The 
rate is always cut off at the dilepton threshold by phase 
space H). 
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FIG. 2. The total integrated rate for dimuons from a 
hadronic gas of T = 50, 100 and 150 MeV. The solid line is 
from the virial expansion and the dotted line from the Haag 
expansion. 

For a parameterization of Fy, we used the common 
Breit-Wigner form with a varying width as found in JIJ] 
and a small correction for large q 2 to fit the data Jig]. 



F v (q 2 ) = 
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with a non-zero width only for q 2 > 4m 2 [ fj"7| 

2 - 4m 2 x 3/2 



TJq 2 ) = 0.149^ 
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and 7 = 0.4. 

For the terms in (||) depending on Im Ily, our results 
will be sensitive to the momentum dependence of the 
width. The KSFR relation implies H 



F v (q') = l + ^U v (q^ + 0[ -L 
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so we have fixed II y from Fy with the l// 4 terms omit- 
ted. Fig. 3 shows this fit against the data taken from the 
recent compilation by Huang p2| . Using the fit instead 
of the data produces a 10% uncertainty at the rho pole 
where the fit is worst. On the log plot of Fig. 1 this is a 
negligible difference. 
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FIG. 3. Data for the vector and axial spectral densities and 
the fits we used. 

The axial channel is dominated by the al(1260). A 
Breit-Wigner fit for FY a works well with (/^^a) = 
(190, 1210) MeV through the substitution m 2 + 717 2 -> 
/|, Trip -» m A , and T p (q 2 ) -> T A (q 2 ) in O. In @J|), the 
threshold is moved to 9m 2 with (m p , 4m^) — * (tua, Qrn 2 ) 
and the overall constant is changed so that TA(m A ) = 0.4 
GeV. The comparison of this fit to the data is also shown 
in Fig. 3. We note that fixed-width calculations tend to 
overestimate the tail of the al in the low mass region by 
a large factor. 

In Fig. 4, we show the non-integrated results for back- 
to-back (q = 0) dielectrons with T = 50, 100, 150 MeV. 
As in the fully integrated rate, our result and the Haag 
expansion are practically identical for the lower tempera- 
tures and, as the temperature increases, the region below 
the two pion threshold increases with respect to the Haag 
rate. 

The shape of the virial expansion curves in Figs 1, 
2, and 4 can be explained as follows. Above the two- 
pion threshold, the first term of eq. (Jfy dominates giving 
the rho peak. As the temperature increases, the Im Ily- 
terms in eq. (Q) start to play a part and, since they give a 
negative contribution to the rate, decrease the rho peak 
by about 10% for T = 150 MeV. The Im U v terms are 
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proportional to T p (q 2 ) and as a consequence all vanish 
below the two-pion threshold. This part is fully consis- 
tent with the Haag rate for all temperatures. 

The main difference comes from the inclusion of the 
Im 11,4 term. It is the only contribution below the two- 
pion threshold. The reason it picks up the al pole contri- 
bution even at low q 2 is because, in contrast to the result 
in the chiral limit, the axial spectral density is integrated 
over all momentum in the thermal averaging. This weak- 
ens the contribution from the (k + q) 2 prefactor in eq. (Q) 
therefore allowing the 1/q 2 dependence in (jl]) to domi- 
nate at low q 2 . Integrating the rate pronounces this even 
more and is only cut off at the dilepton threshold by (0). 
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m,), for T = 100-150 MeV. 

For dileptons of mass above the al pole, our emission 
rate is also larger than the one expected from the Haag 
expansion due to the contribution from the continuum 
in the rho channel as well as the effects from the al. 
At T — 150 MeV the enhancement is about an order of 
magnitude in the dielectron rate. 

If we were to assume an acceptance factor of about 
1/100 (corresponding to the experimental cuts in q± and 
y), a lifetime for the thermal gas of 5 fm/c, and a normal- 
ization of 2n (due to charged particle multiplicities), our 
dielectron rate distribution at a temperature of 100-150 
MeV is similar in shape and magnitude to the physically 
measured dielectron rate by the CERES collaboration [jlj . 
In addition, our dimuon rate distribution for a tempera- 
ture of 50-100 MeV is similar in shape to those measured 
by the HELIOS3 collaboration || prior to the experimen- 
tal acceptance cuts. Of course, these are purely heuristic 
statements that should motivate a more detailed calcu- 
lation using a transport equation. 
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FIG. 4. The back-to-back rate for dielectrons. The virial 
expansion is the solid line and the Haag expansion is the dot- 
ted line. 

This enhancement can be physically understood as fol- 
lows. Since 2 Im W F (g) = W(c?) +W(— q), the low-mass 
region is enhanced by processes such as X — ► ne + e~ 
through W(— q). These are reminiscent of the Dalitz 
decays. This description, however, is very qualitative 
since the relationship between the physical reactions in 
the emission rate and the imaginary part gets blurred by 
the virial expansion. In comparison with most dilepton 
rate calculations [p| ]n^Ji9| , |2l| , our results show a substan- 
tial enhancement in the low mass dilepton region (2-4 
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FIG. 5. The photon emission rate from a hot hadronic gas. 
The solid and dotted lines were calculated with a pion chem- 
ical potential of and 100 MeV respectively. 

Fig. 5 shows the photon emission rate versus the pho- 
ton energy for T = 50, 100, 150 MeV. The solid line is 
the result following from the virial expansion. As we 
have noted, there is no contribution from the Haag ex- 
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pansion to the order quoted. We note that the the bulk 
of the low mass dileptons stem primarily from the pho- 
ton source. Our results are an order of magnitude above 
those obtained using two-body reactions with thermal 
p's and al's |2(],|2l| at the peak of the curve. However, 
above the al mass, our results are in agreement with 
others. We note that the use of a constant width for Ta 
increases the rate by almost an order of magnitude for 
all temperatures. 

Other calculations |2lJ tend to underestimate the pion 
electric polarizability by almost an order of magnitude 
which could point to the source of their result being 
smaller. In our case, the result is constrained by the 
data for 11^. 

The role of the pion chemical potential \i v may be as- 
sessed through the substitution n(uj) = (e"/ T — — > 
(e^-MiO/r — l)" 1 in the above formulas. A non-zero 
chemical potential influences all terms which involve pi- 
ons. This includes the entire Haag rate. In contrast, the 
first term in the virial expansion, which is its dominant 
contribution around the rho pole, is not influenced. 

Fig 6. shows the rates for = 100 MeV which de- 
scribes a dilute system of pions with a mean pion sep- 
aration of about 2 fm at T = 100-150 MeV |IlJ. The 
dilepton integrated rate following from the virial expan- 
sion is enhanced in the low mass region by about a factor 
of five. There is no change near the rho pole and above 
1 GeV it is again slightly enhanced. The entire Haag 
rate, on the other hand, is increased by about a factor 
of 10. This is a major qualitative difference in the two 
analyses. This further enhancement points at the im- 
portance of pions off-equilibrium, thereby motivating a 
more quantitative analysis in the context of kinetic the- 
ory. The photon rates are systematically enhanced as 
shown by the dotted lines in Fig. 5. In the temperature 
range T = 100-150 MeV, the enhancement is about a 
factor of two throughout the photon energy spectrum. 

6. To summarize, we have argued that for a hadronic 
gas at temperatures T < the dilepton and photon 
emission rates were constrained by data from electropro- 
duction, tau decays, radiative pion decay and two-photon 
fusion reactions. Through these constraints, the rates 
were assessed using data for the various vacuum corre- 
lation functions and carried out to the first two terms 
in the virial expansion (@). For completeness, we have 
compared the results following from the Haag expansion. 

Through chiral reduction formulas, our construction 
shows the interdependence between the dilepton and pho- 
ton rates and other low energy processes. We note that 
the chiral reduction formulas follow solely from broken 
chiral symmetry with explicit (2,2) breaking. They are 
a direct generalization of current algebra results from 
threshold to general momenta. In this respect, our anal- 
ysis is different from finite temperature arguments using 
PCAC in the chiral limit jL2|,[l3). Our analysis was con- 
fined to the isovector part of the electromagnetic current, 



and could be improved at higher mass (phi region) by in- 
cluding the isoscalar contribution. 
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FIG. 6. The dielectron integrated rates for T = 50, 100, 150 
MeV and a finite pion chemical potential /i^ = 100 MeV. The 
solid line is from the virial expansion. The dashed line is from 
the Haag expansion. For comparison, the virial expansion 
result at zero chemical potential is shown as the dotted line. 

Since dilepton or photon emission from heavy-ion col- 
lisions involve a dynamical integration over the real-time 
history of the hadronic gas that is expected to expand 
and cool, our results are only suggestive. However, they 
already show a thermal pion gas in equilibrium yields 
dilepton and photon rates that are enhanced in the low 
mass region compared to most of the calculations using 
specific reaction processes. The inclusion of a pion chem- 
ical potential enhances the entire photon rate and the low 
and high mass dilepton rates without affecting the rate 
distribution around the rho mass. These results are rel- 
evant for present dilepton measurements at CERES and 
HELIOS3 as well as future dilepton and photon measure- 
ments at RHIC and LHC. 
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